A NOTE ON ARTIN'S CONSTANT 
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We suggest a representation of Artin's constant, which conjecturally 
describes the density of prime p such that "generic" g G Z is primit- 
ive modulo p. Namely, A = lim^ R k (N) for R k (N) = §£^E$, 
where the summation is over first N prime pi, k G Z + . The clas- 
sical summation formula is as follows: A = \im.N^ 00 'E(N), where 
£(iV) = ^ ^ ^tTi . The changes needed for arbitrary g are addressed 
in Theorem 1, a good exercise in basic analytic and algebraic number 
theory. The same procedure can be applied to other number-theoretic 
constants like A (see, e.g., [Ni]). In Theorem 2, we demonstrate how 
it works for the Stephens constant and for Artin's constants of higher 
ranks (for the density of prime p such that a given set of "generic" 
integers generates Z* ) . 

The following three features of this approach vs. the summation 
formulas are worth noticing. 

1) The restricted summation suggested by P. Moore to make the X— 
formula matching the right heuristic density for arbitrary g G Z gives 
the desired answer in our approach only when g is not a pure (odd) 
power in Z. Otherwise, nontrivial rational multiplicative corrections 
occur; they are calculated in Theorem 1, (ii). 

2) The p k — terms in the denominator and numerator of Rk(N) do not 
influence the limit, which can be heuristically associated with switching 
from primitive roots in Z p to those in (Z/ (p k+1 ))*. The extra p k — factors 
disappear (cancel) in the corresponding summation formula. When 
k — 0, our R— formula for A (without restricting the summation) follows 
from [Pi]. 

3) The R— formulas oscillate significantly around A (and the other 
constants). The magnitude of oscillations increases as k grows; see 
Figure 1. Representing Rk(N) = Yl^Li w i > the weights Wi change 
from 0( l °i+ 2 ) for small i to O( ^ log ^ + ) for i ~ N. Thus, large p 
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receive greater weights in our approach, especially when k is large, 
which increases the range of oscillations. 

1. Brief history. Artin's primitive root conjecture states that given 
an integer g, possibly negative but not a perfect square in Z, the number 
Pjv(g) of prime p among pi = 2,p 2 = 3, . . . ,pn such that g is primitive 
modulo p approaches asymptotically A(g)N as iV — > oo for 

(1) A(g) = A h A d , A d = (l-MMI) II ^Ta II 



p — 'i p 1 — p 

p\d,p\h p\d,p\h 



* = IIP " ^Tf) IE* - J^ny). where 

p|ft ^ p\h ^ ' 

d = Discriminant (Qfy 7 ^]), g = g% for g a E Z and maximal h G N. 

Note that /x( | cZ| ) = and, respectively, A(g) = Ah if and only if the 
discriminant d is not from 1 + 4Z. If h = 1 for such g, then equals 

(2) Artin's constant = A = TT (1 1— — r). 

p(p — 1 

prime p ' 

According to [St], Artin's conjecture was finalized around 1965. In 
1967, it was deduced by Hooley [Ho] from the generalized Riemann hy- 
pothesis for the fields K m = Q[Cm> 9 l ^ m ] f° r squarefree m. See [Mol] for 
a comprehensive introduction (including some recent developments). 
See also [Mu, Le]. 

Artin's heuristic approach to this conjecture was based on the expec- 
tation that events "prime p does not split completely in K q for prime q" 
are independent (subject to later qualitative and quantitative correc- 
tions). For instance, one can expect that A equals limjv^oo VN{g)/N if 
g is generic as far as primitive roots modulo prime p are concerned. It 
leads to the following heuristic summation formula for Artin's constant: 

' 4>(Pi-i) 



(3) A = lim £(iV) , £(A0 = 1 

i=i ^ % 

which can be checked unconditionally, without any reference to Artin's 
conjecture. See [LL], [Mol]. P. Moore extended it to arbitrary A(g) by 
switching to p in this summation such that g is a quadratic nonresidue 
modulo p and (p — 1, h) = 1 (heuristically, it makes sense); see below. 
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As Lehmers wrote, the convergence in (3) is "discouragingly slow" 
(they considered prime numbers p < 1500000). It remains very slow 
when prime numbers in much greater ranges are considered, generally, 
no better than the (conjectural) convergence of Vn(2)/N to A = A(2); 
cf. Table "Artin's constant estimates" from [Si] (p < 10 14 ). 

2. Main Theorem. A refine version of this heuristic approach is from 
[Mol] (for any integer g): 

1 A R H 

( 4 ) *! im Af5>*fo) = ( = j im Vn(9)/N), where 

i=l 



^(p)^2^-^ for (?) = -l and (p-l,h) = l, 
p — 1 \p/ 



tfgip) = otherwise . 
Heuristically, it is equally reasonable to expect that 

E£i ^ - 1) 



(5) A = lim 

N -°° Elife-i) 

for sufficiently general g. Switching here to the restricted summation 
from (4), we come to the following theorem. 

Theorem 1. (i) For arbitrary integers k > and g (possibly, negative), 

(6) A(g) = lim ^% P ^ ~ '\ 

provided that h = 1, where the summation £' is over prime pi such that 
(V) =-1 and ( Pi -l,h) = l. 

(ii) If h > 1, then the limit in the r.h.s. of (6) equals 

(7) A(g) = A h A d , where A h = TT(1 - ) , 

Ad = 1 if d\h and A& = A^ otherwise . 

Sketch of proof . We follow [LL, Mol], restricting ourselves with (5); 
see also [Pi]. Coupling the generalized Landau formula from [SZ] 

Epm ^ (i+faa^ m > o, 

*— ' 1 + m log x 
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with the classical estimate 

ir(x,d,l) = 1 + 0(1/ log a;) 
7r(x) (j)(d) 

where ir(x,d, 1) is the number of prime numbers p < a; in 1 + (i N 
(n(x) = ir(x, 1, 1)), one arrives at: 

7i^(x,d, 1) _ l + 0(x m /logx) 
( ' ttM(x) ~ 0(d) ' 

where 7P m )(x, d, 1) = X]p<zP m over P r i me P £ 1 + dN. Then, 

5>>-i) = £5>-i)^ 

p<x p<x d\p—l 

Finally, 

E*.(p-i) ~* ' 

We will omit the arguments (from the basic algebraic number theory) 
that give the rational corrections for arbitrary g. □ 

In (5,6), only the leading powers of pi matter; for instance, one can 
take here the ratio 0(0(p fc+1 ))/ J2 P 4>{p k+l )i the heuristic probabil- 
ity for g being a primitive root modulo p k+1 over all prime p, which 
leads to the same Artin constant. 



3. Further examples. Let us apply the same procedure to Stephens' 
constant S(a, b) and the higher rank Artin constants. The former is 
defined for a given pair a, b G Q* such that a r b s = 1 r = 0, s = 
for r, s G Z; it describes the density of prime p such that b = a m mod p 
for some m G Z. Modulo the generalized Riemann hypothesis, it equals 

oo 

def TT/1 P* \ 



;il) C ab S for 5 = ] 



where the factors C a b are rational [S, MS]. These factors were calcu- 
lated explicitly in [MS] under the condition that the group Q*/ (a, 6,-1) 
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is torsion free. Assuming that a, b are "random" , the heuristic proba- 
bility P3(p) that b = a m mod p for some m can be readily calculated: 

(12) FS{p) = S{p)/{p-l) 2 , where S(p) = ^ 

d\ P -i 

= J^J — for prime factorization p — 1 = J^J qy> . 

j=i qj + i=i 

We naturally omit the prime numbers that divide the numerators or 
denominators of a, b. 

The rank r Artin constant A r (gi, . . . , g r ) describes the heuristic den- 
sity of prime p such that a given set of nonzero integers {gi, . . . ,g r } 
(or rationals) generates Z* . See [CP] . Its "generic" value (modulo the 
generalized Riemann hypothesis) is as follows: 

OO i 

def 



as) A r =n i 



p[ (Pi 



The corresponding probability at p for "random" gi, . . . ,g r equals 

m 

(14) PA r {p) = A r (p)/(p - IT for Mp) = \{{qf - qfi^) 

in terms of the prime factorization p — 1 = rij=i Qj 3 • 

The summation (unconditional) limiting formulas £1X6 cLS follows: 

iV AT 

(15) S = lim ^^iSfe), A- = lim ^PA r fe). 

N^oo iv ^ N^oo iv ^ 

i=l j=l 

The proof of the following theorem is similar to that of Theorem 1. 
Theorem 2. For an arbitrary integer k > 0, 

(16) 5 = lim ^/=^ S ^ , 

(17) A r = lim ^iP?A(Pi) 

where the summation is over consecutive prime numbers pi. □ 
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4. Numerical aspects. For the constants C considered above, we 
plot ~^p- — 1, shown blue-thin, for the classical summation £(iV) and 

— 1, which are red-thick, for our ratio approximations R(N); for 
instance, £(AQ is from (3) for Artin's A. The range is iV < 1000M. 

DEVIATION FROM A, STEP = 5000000 PRIMES 
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DEVIATION FROM A fork=l, STEP = 5000000 PRIMES 



0.00001 - 



5.x 10 



-5.x 10" 
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Figure 1. Deviation from A as k = 0, 1 for 1000M primes. 
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Figure 1 compares the stabilization of (3), blue-thin, to the Artin 
constant and the stabilization of (5), red-thick, as k = 0, 1. The func- 
tion S(A r ) — A remains positive in the range N < 1000M; R^(N) os- 
cillate around Artin's constant A » 0.37395581361920228805 (the zero 
level of this graph). The amplitude of oscillations become larger for 
Rk = i vs. Rk=0i however the graphs are very much similar. 

We note that the best way to calculate A and similar constants is 
based on the known product formulas in terms of ((n) for integers 
n > 1. 




Figure 2 shows the convergence of (4) and (6) to A(g — 5) for k = 0; 
notice that the blue-thin curve remains beyond the red-thick one in 
this range. 

The last two plots show the graphs for the Stephen constant S ~ 
0.57595996889294543964, Figure 3, and the rank 3 Artin constant A 3 » 
0.85654044485354217443, Figure 4; here k = 0. The convergence rate 
and other features of these four graphs are similar to those for A and 
A(5). There is striking (qualitative) similarity of these two figures, 
including the oscillations, although the convergence rate in Figure 4 
(for A 3 ) is significantly (almost 10 times) greater than in Figure 3. 
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DEVIATION FROM S, STEP = 5000000 PRIMES 
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Figure 3. Deviation from S for 1000M primes. 
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DEVIATION FROM A[r=3], STEP = 5000000 PRIMES 
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Figure 4. Deviation from A 3 (rank=3) for 1000M primes. 
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Qualitatively, the behavior of R{N) for large N can be evaluated 
following [Pi]. For instance, the functions \R(N)/A — 1| in Figure 1 
must be no greater than C m (log N)~ m for any fixed m > and proper 
constant C m (depending on k) as N » 0. Generally speaking, C m 
can be estimated in terms of (the order of) N, but we will not discuss 
it. Indeed, the graphs of R(N)/A — 1, red-thick from Figure 1, look 
like 0((logiV)~ 4 ) in the range N < 1000M. 

The nature of oscillations of the functions R(N) around the corre- 
sponding constants remains unclear. 

Acknowledgements. I am grateful to Zeev Rudnik for valuable com- 
ments. I am very thankful to Pieter Moree for reading the note, sug- 
gesting interesting questions toward comparing our approach with the 
results from [Mo2, Mo3] and for the reference to [Pi]. 
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